Optical determination of the relation between the electron-boson coupling function 
and the critical temperature in high cuprates. 
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We take advantage of the connection between the free carrier optical conductivity and the glue 
function in the normal state, to reconstruct from the infrared optical conductivity the glue-spectrum 
of ten different high-Tc cuprates revealing a robust peak in the 50-60 meV range and a broad con- 
tinuum at higher energies for all measured charge carrier concentrations and temperatures up to 290 
K. We observe that the strong coupling formalism accounts fully for the known strong temperature 
dependence of the optical spectra of the high Tc cuprates, except for strongly underdoped samples. 
We observe a correlation between the doping trend of the experimental glue spectra and the critical 
temperature. The data obtained on the overdoped side of the phase diagram conclusively excludes 
the electron-phonon coupling as the main source of superconducting pairing. 



I. INTRODUCTION. 

The theoretical approaches to the high Tc pairing 
mechanism in the cuprates are divided in two main 
groups: According to the first electrons form pairs due 
to a retarded attractive interaction mediated by virtual 
bosonic excitations in the solid^-^''^''*'^. These bosons can 
be lattice vibrations, fluctuations of spin-polarization, 
electric polarization or charge density. The second group 
of theories concentrates on a pairing-mechanism entirely 
due to the non-retarded Coulomb interaction^ or so- 
called Mottnessi. Indeed, optical experiments have found 
indications for mixing of high and low energy degrees of 
freedom when the sample enters into the superconducting 
stateMOOdi. 

An indication that both mechanisms are present was 
obtained by Maier, Poilblanc and Scalapino^, who 
showed that the 'anomalous' self-energy associated with 
the pairing has a small but finite contribution extend- 
ing to an energy as high as U , demonstrating that the 
pairing-interaction is, in part, non-retarded. The exper- 
imental search for a pairing glue will play an essential 
role in determining the origin of the pairing interaction. 
Aforementioned glue is expressed as a spectral density 
of these bosons, indicated as a^F((jj) for phonons and 
for spin fluctuations, here represented as the gen- 
eral, dimensionless function n(a;). An important conse- 
quence of the electron-boson coupling is, that the energy 
of the quasi-particles relative to the Fermi level, ^, is 
renormalized, and their lifetime becomes limited by in- 
elastic decay processes involving the emission of bosons. 
The corresponding energy shift and the inverse lifetime, 
i.e. the real and imaginary parts of the self-energy, are 
expressed as the convolution of the 'glue- function' n(ijj) 
with a kernel K{^,uj,T) describing the thermal excita- 
tions of the glue and the electrons^ 



In the absence of a glue and of scattering off impuri- 
ties the effect of applying an AC electric field to the 
electron gas is to induce a purely reactive current re- 
sponse, characterized by the imaginary optical conduc- 
tivity Ana {to) = iUp/uj, where the plasma frequency, Wp, 
is given by the (partial) f-sum rule for the conduction 
electrons. The effect of coupling the electrons to bosonic 
excitations is revealed by a finite, frequency dependent 
dissipation, which can be understood as arising from pro- 
cesses whereby a photon is absorbed by the simultaneous 
creation of an electron-hole pair and a boson. As a result, 
the expression for the optical conductivity in the normal 
state, 
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now contains a memory function or optical self 
energyi^ii^. A particularly useful aspect of this repre- 
sentation is that M{uj) follows in a straightforward way 
from the experimental optical conductivity. The optical 
self-energy is related to the single particle self-energies 
by the expression^'' 
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The central assumption in the above is the validity of 
the Landau Fermi-liquid picture for the normal state. 
The aforementioned strong coupling analysis is there- 
fore expected to work best on the overdoped side of the 
cuprate phase diagram, where the state of matter appears 
to become increasingly Fermi liquid like. If antiferromag- 
netism is necessary to obtain the insulating state in the 
undoped parent compounds, as has been argued based 
on the doping trends of the Drude spectral weightJi, the 
strong coupling analysis may in principle be relevant for 
the entire doping range studied. However, in the limit of 
strong interactions aforementioned formalism needs to be 
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extended, e.g. with vertex corrections, and it eventually 
breaks down. We therefore define the function n(w) as 
the effective quantity which, in combination with Eqs. [T] 
and[31 returns the exact value of M{uj) for each frequency. 
Defined in this way n(a;) captures all correlation effects 
regardless whether the system is a Fermi- liquid or not. 
This becomes increasingly relevant when the doping is 
lowered below optimal doping. 

Here we take advantage of the connection between the 
temperature and frequency dependent conductivity in 
the normal state and the glue-spectrum to test experi- 
mentally the consequences of the standard approach, to 
check the internal consistency of it, and to determine 
the range of doping where internal consistency is ob- 
tained. For a d-wave superconductor, the momentum 
dependence is essential to understand the details of the 
pairing. This, of course, is difficult to handle for opti- 
cal spectroscopy which is inherently a momentum inte- 
grated probe. Nevertheless, optical spectra provide the 
important information on the energy scale of the bosons 
involved and on the doping and temperature evolution. 
The paper is organized as follows. In section |TT] we show 
that the temperature dependence of the optical spectra of 
the cuprates is well described within the strong coupling 
formalism described above. In section [1111 we present the 
n(ti>) functions for 10 different cuprates. These n(w) 
functions are used in section IIVI to estimate critical tem- 
peratures and section |V] discusses the implications of 
these results with regard to the pairing mechanism in 
the cuprates. Finally, in section IVII we summarize our 
results. 



II. INTERNAL CONSISTENCY CHECK OF 
THE STRONG COUPLING FORMALISM. 

In order to test whether the strong coupling analysis 
is applicable to the cuprates we start with an important 
test of its internal consistency: (i) we invert the data 
at 290 K to obtain n(w), (ii) we use this n(u;, 29QK) to 
predict the optical spectra at lower temperatures. If the 
prediction faithfully reproduces the experimental spectra 
at these temperatures, we have a strong indication that 
the electronic structure and its evolution as a function 
of temperature are to a good approximation within the 
realm of strong coupling theory. We use a standard least 
squares routine to fit a histogram representation of n(a;) 
to our experimental infrared spectra (see Appendix). The 
quantity IV{lo) is shown in Fig. [1] for optimally doped 
HgBa2Cu04+5 (Hg-1201)is for T = 290 K, together with 
the optical self energies calculated from this function at 
three different temperatures. For 290 K the theoretical 
curve runs through the data points, reflecting the full 
convergence of the numerical fitting routine. 

It is interesting to notice, that the shoulder at 80 meV 
in the 100 K experimental data is reproduced by the same 
n(w) function as the one used to fit the 290 K data. It 
can be excluded that this shoulder is due to the pseudo- 
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FIG. 1: Experimental optical self energy of HgBa2Cu04+i 
for 3 selected temperatures (open circles). The solid curve 
at 290 Kelvin is obtained from a fit of n(a;) , shown as the 
dashed surface. The solid curves at 100 K and 20 K were 
calculated with the same n(w) function corresponding to 290 
Kelvin. This proves that the self energy feature between 80 
and 100 meV (a shoulder at 100 K and a peak at 20 K) is 
caused by the prominent peak in n(cij) at approximately 60 
meV. The sharpening of this feature at low temperature is due 
to the superconducting gap, an aspect not captured by Eq. [3] 
and therefore not reproduced in the calculated solid curves. 
In the inset the gap-induced sharpening is illustrated by the 
optical self energy without (black) and with (red) a 15 meV 
superconducting gap, calculated using Allen's relatiorji^. 



gap, since a gap is certainly absent for temperatures as 
high as 290 K. The shoulder is therefore entirely due to 
coupling of the electrons to a mode at approximately 60 
meV. On the other hand, the considerable sharpening 
of this feature for temperatures lower than 100 K finds 
a natural explanation in the opening of a gap, as illus- 
trated in the inset of Fig. [1] We see that for 100 K the 
theoretical prediction also runs through the experimental 
data points. In other words, the strong temperature de- 
pendence of the experimental optical spectra is entirely 
due to the Fermi and Bose factors of Eqs. [T]and[3l 

This example confirms the close correspondence be- 
tween the features in M(uj) and in n(a;) pointed out in 
Ref. [l^. In particular the broad maximum in M(a;) has 
its counterpart in the high intensity region of n(u;) ter- 
minating at 290 meV. The internal consistency is there- 
fore demonstrated by the fact that the large tempera- 
ture dependence of the optical spectra is fully explained 
by the strong coupling formalism. This consistency was 
obtained for all samples, except for the most strongly 
underdoped single layer Bi2201 sample. 



III. ELECTRON BOSON COUPLING 
FUNCTION. 

As summarized in Fig. [21 we have analyzed previously 
published optical spectra of 10 different samples belong- 
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FIG. 2: Electron-boson coupling function Tl{uj) for Bi-2201 at four different charge carrier concentrations (10 K, 100 K, 290 K), 
Bi-2212 at four charge carrier concentrations, and optimally doped Bi-2223 and Hg-1201 (100 K, 200 K, 290 K). The samples 
are ordered from underdoped to overdoped (left to right) and low to high Tc (top to bottom) 



ing to different families of materials, i.e. optimally doped 
Hg-120li^ and Bi2Sr2Ca2Cu30io+5 (Bi-2223)ii, as well 
as four Bi2Sr2CaCu208+5 (Bi-2212) crystals^ with dif- 
ferent hole concentrations. In addition, we analyzed data 
for four Bi2Sr2Cu2 06+5 (Bi-2201) crystals with different 
hole concentrations2£.. 

Excellent fits were obtained for all temperatures, but 
the n((jj) spectra exhibit a significant temperature depen- 
dence, in particular at the low frequency side of the n(cj) 
spectrum. Since all thermal factors contained in Eqs. 
[T] and [3] are, in principle, folded out by our procedure, 
the remaining temperature dependence of n(aj) reflects 
the thermal properties of the 'glue-function' itself. Such 
temperature dependence is a direct consequence of the 
peculiar DC and far infrared conductivity, in particular 
the T-linear DC resistivity and uj/T scaling of T<t{llj,T) 
at optimal doping^ ^. For the highest doping levels both 
n(ci;) and its temperature dependence diminish, which is 
an indication that a Fermi liquid regime is approached. 
The most strongly underdoped sample, Bi-2201-UDO, ex- 
hibits an upturn of the imaginary part of the experimen- 
tal optical self-energy for w ^ 0. This aspect of the data 
can not be reproduced by the strong coupling expression. 



resulting in an artificial and unphysical peak at w « of 
the fitted II(cj) function. 

We observe two main features in the glue- function: A 
robust peak at 50-60 meV and a broad continuum. The 
upper limit of n(w) is situated around approximately 300 
meV for optimally doped single layer Hgl201, and for the 
bilayer and trilayer samples. The continuum extends to 
the highest energies (550 meV for the single-layer samples 
and 400 meV for the bilayer) for the weakly overdoped 
samples, whereas the continuum of the strongly doped 
bilayer sample extends to only 300 meV. There is also a 
clear trend of a contraction of the continuum to lower en- 
ergies when the carrier concentration is reduced. Hence, 
part of the glue function has an energy well above the 
upper limit of the phonon frequencies in the cuprates (~ 
100 meV). Consequently the high energy part of tl{uj) re- 
flects in one way or another the strong coupling between 
the electrons themselves. 

The most prominent feature, present in all spectra re- 
produced in Fig. [21 is a peak corresponding to an av- 
erage frequency of 60 ± 3 meV at room temperature 
(see Appendix for an estimate of the error bar) . Perhaps 
the most striking aspect of this peak is the fact that its 
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TABLE I: Strong coupling parameters of the ten compounds. 
The hole-doping is indicated on the first line. From left to 
right: Bi2201 (columns 1 to 4), Bi2212 (columns 5 to 8), 
Bi2223 (9th column) and Hgl201 (column 10). On rows 6-9 
we indicate the partial coupling constants Tc's obtained when 
the Tl{uj) spectra are separated in a contribution from the 
peak (pk, lo < 100 meV) and from the continuum (cnt, cu > 
100) meV. All values are listed for room temperature. 

energy is practically independent of temperature (up to 
room temperature) and sample composition. Moreover, 
the intensity and width are essentially temperature inde- 
pendent. While our results confirm by and large the ob- 
servations of Hwang et al. in the pseudo-gap phase^^^^i, 
the persistence of the 50-60 meV peak to room tempera- 
ture has not been reported before for these compounds. 
However, Collins et al obtained excellent fits to their in- 
frared data of YBa2Cu307 at 100 K and 250K using for 
both temperatures the same a^F^uj) spectrum with a 
peak at ~ 35 meV and a continuum extending up to 
300 meV. The 50-60 meV peak which we observe, arises 
most likely from the same boson that is responsible for 
the 'kink' seen in angle resolved photoemission (ARPES) 
experiments along the nodal direction in k-space at ap- 
proximately the same energy^^i^i^. The peak-dip-hump 
structure in the tunneling spectra (STS}^>^'^" has also 
been reported at approximately the same energy. 



the pairing interaction is obtained by taking the average 
of ln(a;) weighted by electron-boson coupling functional, 

/•OO 

ln(c:)) = 2A"W uj-^Ild{oj)ln{uj)duj. (6) 
Jo 

To apply Eq. 0] to our experimentally measured TI{lu) , 
we would need to map this function on the d-wave pair- 
ing channel. Boson fluctuations below a certain critical 
frequency act as pair breakers, as has been shown by Mil- 
lis, Varma and Sachdev2? in the case of spin-fluctuation- 
mediated d-wave superconductivity. Clearly, it is not 
possible to separate pair-breaking from pair-forming con- 
tributions to H(lj) in an unambiguous way. To proceed 
we assume that the full n(a;) function contributes favor- 
ably to the pairing. This means that our results overesti- 
mate the critical temperature. In Table [1] we indicate the 
total coupling constant. A, and logarithmic frequency, (D, 
for the room temperature n(a;) spectra. The coupling 
strength shows a strong and systematic increase with de- 
creasing hole concentration, which probably requires a 
theoretical treatment beyond the strong coupling expan- 
sion. At the same time we see that J) shows the opposite 
trend. 

An estimate of Tc, using the experimental values in- 
dicated in Table IH gives values in the 100-200 K range. 
The critical temperature can also be calculated straight- 
forwardly from the s-wave Eliashberg equations^ when 
Il{uj) is known. As shown in Fig. [3l the T^s are in 
the 150-300 K range, and they correlate with the experi- 
mentally observed doping trends of Tc- The dome-shaped 
trend in the calculation is a consequence of the increasing 
energy scale of H(ix') and the decreasing overall coupling 
constant as a function of doping. 

V. IMPLICATIONS FOR THE PAIRING 
MECHANISM. 



IV. CRITICAL TEMPERATURES 

One of the most important issues in the field of high Tc 
is the question whether pairing is caused by the exchange 
of virtual bosons. These processes are described by a 
bosonic density of states function closely related to n(a;) . 
If the electron-electron interaction occurs uniquely in the 
d-wave channel, the superconducting critical temperature 
follows from the usual relation 

Tc-0.83wexp(-(l + A<i)/Arf), (4) 

where A^ is the coupling constant in the d-wave pairing 
channel is 

/"OO 

Ad = 2 / Ud{uj)/ujduj, (5) 
Jo 

and Ild{uj) is the d-wave electron-boson coupling func- 
tion. The effective frequency of the bosons responsible for 



We take this analysis a step further by calculating Tc 
from the glue spectra below 100 meV (Hp^ ) and above 
100 meV (Hcnt). The resulting coupling constants and 
Tc's are indicated in Table [J On the underdoped side 
Ilcnt vanishes, and Tc is given only by the coupling to the 
intense 50-60 meV peak in n(w), but in this doping range 
we have to be careful with the interpretation of our re- 
sults. As mentioned in the introduction our H{i-u) spectra 
represent effective coupling functions, which may contain 
effects arising from features not captured by the strong 
coupling equations [T] and [3) If, for example, a pseudogap 
opens in the electronic spectrum this will affect the shape 
of n(w). These effects likely play a role for the under- 
doped samples, but are not expected to affect much the 
room temperature values, indicated in Tablc|T]and Fig. [31 
On the contrary, the larger temperature dependence seen 
for underdoped samples in Fig. [2] may well be a result of 
the opening of a pseudogap. For the overdoped samples 
the Tc's calculated from Hp^ are smaller than the exper- 
imental values. For example, for Bi2212 with the highest 
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FIG. 3: Experimental critical temperature (open symbols) 
and Tc's calculated in the Eliashberg formalism (closed sym- 
bols) using the experimentally measured Ti{ui) of Fig. [2] at 
290 Kelvin as input parameters. 



doping Upk gives only Tc < 20K, whereas Ilcnt gives 160 
K, implying that the glue-function above 100 meV is of 
crucial importance for the pairing- mechanism. Since only 
electronic modes can have such high energies, an impor- 
tant contribution to the high Tc mechanism comes ap- 
parently from coupling to electronic degrees of freedom, 
i.e. spinii^iii^iii^ or orbital current fluctuations^. 



VI. SUMMARY 

In summary, the n(w) spectrum obtained from the op- 
tical spectra of 10 different compounds using a strong 
coupling analysis, is observed to consist of two features: 
(i) a robust peak in the range of 50 to 60 meV and (ii) a 
doping dependent continuum extending to 0.3 eV for the 
samples with the highest Tc. We perform an important 
test of the internal consistency of the strong coupling for- 
malism by showing that the temperature dependence of 
the optical spectra is determined by Fermi and Bose fac- 
tors in the strong coupling expressions. The remaining 
temperature dependence of TI{lo) can therefore be taken 
to indicate that part of the Il{uj) spectrum is electronic in 
origin. We observe an intriguing correlation between the 
doping trend of the experimental glue spectra and the 
critical temperature. Finally we obtain an upper limit 
to the contribution of electron-phonon coupling to the 
pairing of the overdoped samples, which is too small to 
account for the observed critical temperature. 
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VIII. APPENDIX 
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FIG. 4: Experimental reflectivity (red and blue lines) and fit 
curves (black lines) for selected samples and temperatures, 
(a): underdoped non-superconducting Bi-2201 (UDO). (b): 
optimally doped Hg-1201 with n ~ 97 K (OpD97) [Q (c): 
overdoped non-superconducting Bi-2201 (ODO). Weak sharp 
peaks, particularly visible for the strongly underdoped sample 
in panel (a) are due to transverse optical phonons, which we 
do not intend to fit. 

The inversion of Eq's 1 and 2 allows to extract 11 (w) 
from experimental data of the optical conductivity, or 
related optical spectra. The accuracy of the resulting 
n(Lj) spectrum is in practice limited by the convolution 
with thermal factors expressed by Eq's 1 and 2^. Mi- 
croscopic models giving roughly the same n(w) spectra, 
which differ however in the details of the frequency de- 
pendence of this quantity, may therefor provide fits to the 
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directly measured optical quantities, such as infrared re- 
flectance spectra, which at first glance look satisfactory, 
but the remaining discrepancies with the experimental 
spectra may nevertheless be of significant importance for 
the proper understanding of the optical data. It is there- 
fore of crucial importance to test the 'robustness' of each 
fit with regard to the spectral shape of the II (w) function 
imposed by such models. This robustness can be tested 
by including in the fit-routine one or several 'oscillators' 
superimposed on the model function. When the model 
glue function provides a complete description of the elec- 
tronic structure, adding extra oscillators will not result in 
an improvement of the quality of the fit. We have used 
this approach to test functional forms commonly used 
in the literature, in particular the marginal Fermi liq- 
uid (MFL) model^ and the Milhs-Monien-Pines (MMP) 
representation of the spin fluctuation spectrun*^. 

We found that neither of these functional forms de- 
scribe completely the experimental data. In search of 
a more flexible form of n(cj) we used a superposition 
of lorentzian oscillators and found that it could be used 
to describe all available experimental data in a consis- 
tent manner. The resulting n(cj) functions and trends 
are equivalent to those in Fig. [51 From these initial 
tests we concluded that due to the thermal smearing ex- 
pressed by Eq's 1 and 2 our Yl{uj) spectra can only be 
determined with limited resolution. This lead us to the 
use of a histogram representation, where each block in 
the histogram represents a likelihood to find coupling to 
a mode with a well determined coupling strength. For 
the lowest frequency interval {0 < lu < uji) a triangu- 
lar shape was used instead of a block, which is necessary 
to avoid problems with the convergence of the integral 
A = 2 u!^^Il{LL!)duj. In practice the output generated 
by the fitting routine has low intensity in this first inter- 
val, and the triangles are therefore difficult to distinguish 
in Fig. H 

To give an example: The block centered at 55 meV 
seen in the Hg-1201 sample in Fig. [2] has A 1 and 
a width of about 30 meV. Our histogram representation 
implies the presence of a coupling to one or several modes 
between 45 meV and 75 meV with an integrated coupling 
strength of 1. The histograms thus constitute the most 
detailed representation of n(a;) given the precision of our 
experimental reflectivity and ellipsometry spectra. 

Examples of experimental reflectivity data together 
with the fits are shown in Fig. 3] for a selection of repre- 
sentative data sets spanning the entire doping and tem- 
perature range. As the fitted curves are within the limits 
of the experimental noise, further reduction of x^, while 
in principle possible by fitting the statistical noise of the 
data, can not improve the accuracy of the n((jj) functions. 

Starting from a n(u;) function we can calculate the op- 
tical conductivity, which in turn is fed into standard Fres- 
nel expressions to calculate the experimentally measured 
quantities, i.e. reflectivity and ellipsometric parameters. 
The fitting routine is based on the Levenberg-Marquardt 
algorithm and uses analytical expressions for the partial 



derivatives of the reflectivity coefficient R, and the ellip- 
sometric parameters ijj and A relative to the parameters 
describing the fl{to) function. The algorithm is based on 
minimizing a functional which is given by, 

^2 ^ ^ - f{uJi,pi, ...,p„) \ ^ 

z— 1 ^ ' 

where R{uJi) is an experimentally measured datapoint, 
f{LUi,pi, ...,Pn) is the calculated value in this point based 
on parameters pi, ...,Pn and the difference between these 
two is weighed by the errorbar ai determined for R[ijJi). 
For a given set of reflectivity and ellipsometry data at 
one particular temperature, using a standard PC, the it- 
eration takes about 3 hours until convergence is reached. 
The Levenberg-Marquardt least squares method is an ex- 
tremely powerful method to find the minimum of in 
a multidimensional parameter space. To ensure that 
has converged to the global minimum in parameter space 
several tests have been performed, for each individual 
sample and temperature displayed in Fig. where in 
each test the optimization process was started from a 
different set of starting parameters. To give some idea 
of the robustness of our method we will here discuss one 
representative example: optimally doped IIg-1201. 

The models are evaluated based on the minimum found 
for x^. A comparison of Fig. [5] (a-d) shows that the 
MMP model describes better the optical data then the 
MFL model but that they give similar results if we add 
an extra oscillator to these models. Panels [5] (e-f ) show 
the model independent results mentioned above and are 
very similar to the modified MMP and MFL model. The 
models in these last two panels have the same and 
the comparison in Fig. [5^ shows that the histogram rep- 
resentation realistically expresses the uncertainty in the 
position of the low energy peak, while the correspondence 
between the features in both models remains excellent. 
It is interesting that the model with two oscillators is de- 
scribed by 6 parameters, while the histogram represen- 
tation uses 12 parameters. The fact that the fit-routine 
adjusts the latter 12 parameters in such a manner as to 
reproduce the two oscillators, proves that the features 
represented in the righthand panel of Fig. [5] are realistic. 

The models presented in figure [5] allow us to make an 
estimate of the uncertainty in the determination of the 
frequency of the low energy peak. We define the first 
moment of this peak as, 

(uj) = J ujn{uj)duj / J U{uj)duj. (8) 

The variance of (uj) is defined as, 

i=l 

with {uj) the mean of the moments of the spectra pre- 
sented in figure ini and i = 1...N runs over the number 
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FIG. 5: Comparison of several models. The number of data points over which Xj is summed (see Eq. 1) is A*' = 400. The 
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(MMP). In contrast, in the oscillator and the histogram model the was found to be independent of temperature. 
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of spectra used {N — 6 for each temperature). For the 
Hgl201 room temperature spectra presented in figureEK- 
f we find {to) w 60 meV and a 3 meV. This value is 
approximately the same for all samples. In figure [6] we 
present the temperature dependence of the first moment 
of the glue functions presented in figure [H 



FIG. 6: Temperature dependence of (uj) for the Tl{u;) spectra 
presented in figure (5] The error bars are defined in Eq. [51 



^ D.J. Scalapino, E. Loh, J.E. Hirsch, Phys. Rev. B 34, 8190 ^ CM. Varma, P.B. Littlewood, S. Schmitt-Rink, E. Abra- 
(1986). 



8 



hams, A.E. Ruckenstein, Phys. Rev. Lett. 63, 1996 (1989). 
^ A.J. Millis, H. Monien, D. Pines, Phys. Rev. B 42, 167 
(1990). 

* S.V. Shulga, O.V. Dolgov and E.G. Maksimov, Physica C 
178, 266 (1991). 

^ Ar. Abanov, A.V. Chubukov, J. Schmalian, J. Elec. Spec. 

Rel. Phen. 117, pl29 (2000). 
^ P. W. Anderson, Science 316, 1705 (2007). 

P. Phillips Ann. Phys. 321, 1634 (2006). 

* D. N. Basov,S. I. Woods, A. S. Katz, E. J. Singley, R. C. 
Dynes, M. Xu, D. G. Hinks, C. C. Homes, and M. Strongin, 
Science 283, 49-51 (1999). 

® H.J. A. Molegraaf, C. Presura, D. van der Marel, P.H. Kes, 
M. Li, Sctence 295, 2239 (2002). 

F. Carbone, A. B. Kuzmenko, H. J. A. Molegraaf, E. van 
Heumen, V. Lukovac, F. Marsiglio, D. van der Marel, K. 
Haule, G. Kotliar, H. Berger, S. Courjault, P. H. Kes, and 
M. Li Phys. Rev. B 74, 064510 (2006). 
F. Carbone, A. B. Kuzmenko, H. J. A. Molegraaf, E. van 
Heumen, E. Giannini, and D. van der Marel, Phys. Rev. B 
74, 024502 (2006). 

T. A. Maier, D. Poilblanc, and D. J. Scalapino, Phys. Rev. 

Lett. 100, 237001 (2008). 
- m,..T) = J [ + de where n{.) 

and /(e) are the Bose and Fermi-Dirac distribution func- 
tions respectively^^. 
" W. Goetze, P. Woelfle, Phys. Rev. B 6, 1226 (1972). 

An alternative frequently used notation for this memory 

functions^ is -2E°''(a;) = M{u). 
^'^ P.B. Allen, Phys. Rev. B 3, 305, (1971). 
^"^ A. Comanac, L. de' Medici, M. Capone, A.J. Millis, Nature 

Phys. 4, 287, (2008). 

E. van Heumen, R. Lortz, A.B. Kuzmenko, F. Carbone, 

D. van der Marel, X. Zhao, G. Yu, Y. Cho, N. Barisic, M. 
Greven, C.C. Homes, and S. V. Dordevic Phys. Rev. B 75, 
054522 (2007). 

M.R. Norman, A.V. Chubukov, Phys. Rev. B 73, 
140501 (R) (2006). 

E. van Heumen, A.B. Kuzmenko, D. van der Marel, H. 
Eisaki, and W. Meevasana, unpublished. 



■^^ D. van der Marel, H. J. A. Molegraaf, J. Zaanen, Z. Nussi- 
nov, F. Carbone, A. Damascelli, H. Eisaki, M. Greven, P. 

H. Kes, and M. Li, Nature 425, 271 (2003). 

22 J.L. Tallon, J. W. Loram, G. V. M. Wilhams, J. R. Cooper, 

I. R. Fisher, J. D. Johnson, M. P. Staines, and C. Bernhard, 
phys. Stat. sol. (b) 215, 531 (1999). 

23 J. Hwang, T. Timusk, E. Schachinger, J.P. Carbotte, Phys. 
Rev. B 75, 144508 (2007). 

2* J. Hwang, E.J. Nicol, T. Timusk, A. Knigavko, J.P. Car- 
botte, Phys. Rev. Lett. 98, 207002 (2007). 
2^ P.V. Bogdanov, A. Lanzara, S. A. Kellar, X. J. Zhou, E. 

D. Lu, W. J. Zheng, G. Gu, J. -I. Shimoyama, K. Kishio, 

H. Ikeda, R. Yoshizaki, Z. Hussain, and Z. X. Shen, Phys. 

Rev. Lett. 85, 2581 (2000). 
2'^ A. Lanzara, P. V. Bogdanov, X. J. Zhou, S. A. Kellar, D. 

L. Feng, E. D. Lu, T. Yoshida, H. Eisaki, A. Fujimori, K. 

Kishio, J.-I. Shimoyama, T. Noda, S. Uchida, Z. Hussain, 

Z.-X. Shen, Nature 412, 510 (2001). 
2^ W. Meevasana ,N. J. C. Ingle, D. H. Lu, J. R. Shi, F. 

Baumberger, K. M. Shen, W. S. Lee, T. Cuk, H. Eisaki, 

T. P. Devereaux, N. Nagaosa, J. Zaanen, and Z.-X. Shen, 

Phys. Rev. Lett. 96, 157003 (2006). 
28 J. Lee, K. Fujita, K. McElroy, J. A. Slezak, M. Wang, Y. 

Aiura, H. Bando, M. Ishikado, T. Masui, J.-X. Zhu, A. V. 

Balatsky, H. Eisaki, S. Uchida, and J. C. Davis, Nature 

442, 546 (2006). 
2^ G. Levy deCastro, C. Berthod, A. Piriou, E. Giannini, and 

O. Fischer, Phys. Rev. Lett. 101, 267004 (2008) 
3" J.F. Zasadzinski, L. Ozyuzer, N. Miyakawa, K. E. Gray, D. 

G. Hinks, and C. Kendziora, Phys. Rev. Lett. 87, 067005 

(2001). 

3^ B. Mitrovic, H.G. Zarate and J.P. Carbotte, Phys. Rev. B 
29, 184 (1984) 

32 A.J. Millis, S. Sachdev, and CM. Varma, Phys. Rev. B 
37, 4975 (1988). 

33 C.S. Owen and D.J. Scalapino, Physica 55, 691 (1971). 

3" S.V. Dordevic, C. C. Homes, J. J. Tu, T. VaUa, M. Stron- 
gin, P. D. Johnson, and G. D. Gu, D. N. Basov Phys. Rev. 
B 71, 104529 (2005). 



